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Abstract 

The main features of QCD, e.g. confinement, chiral symmetry 
breaking, Regge trajectories are naturally and economically explained 
in the framework of the Field Correlator Method (FCM). The same 
method correctly predicts the spectrum of hybrids and glueballs. When 
applied to DIS and high-energy scattering it leads to the important 
role of higher Fock components in the Fock tower of moving hadron, 
containing primarily gluonic excitations. 



1 Introduction 



The QCD is the (only) internally selfconsistent theory, denned by the only 
scale (parameter the string tension a or Aq C d can be used for this pur- 
pose), which eventually explains fundamental structure of baryons and other 
hadrons and by this some 98% of all the mass in our world. The features of 
QCD are unique in their complexity: confinement, chiral symmetry break- 
ing, string structure of hadrons demonstrating nonperturbative (NP) inter- 
actions, on one hand and on another hand the applicability of perturbation 
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theory for large Q and large momentum processes, where NP effects enter 
as corrections. The most popular theoretical approaches to these phenomena 
look fragmentary, i.e. magnetic monopoles are used for confinement, instan- 
tons for chiral symmetry breaking and pure perturbative expansions for high 
energy processes. 

The situation has improved with the introduction of the QCD sum rules 
PP, where NP contributions are encoded in the form of local condensates. In 
this talk I shall describe a general method which allows to consider all NP 
effects on one ground - with the help of nonlocal Field Correlators (FC) [2] • 
It will be argued that the simplest of these correlators - the quadratic one 
- is dominant and is responsible for confinement, chiral symmetry breaking, 
and the structure of meson and baryon spectra. 

When supplemented with the Background Perturbation Theory (BPT) 
the method allows to treat valence gluons in the confining background [3]. 
This yields the spectrum of hybrids and glueballs in good agreement with 
lattice simulations, and the new perturbation series for high energy processes, 
without IR renormalons and Landau ghost poles. As the new and unexpected 
element, it will be argued that hybrids play a more fundamental role in DIS 
and high-energy scattering - as the building blocks of the colliding hadron 
wave functions. 

The talk is organized as follows. In section 2 FC are introduced and 
confinement is related to their properties. In section 3 the Hamiltonian for 
valence components is written down and properties of meson spectrum are 
discussed. In section 4 this Hamiltonian is extended to the systems with 
valence gluons and spectrum of hybrids and glueballs is discussed. In sec- 
tion 5 the Fock towers are introduced for hadrons and the general matrix 
Hamiltonians is written down. In section 6 the role of higher (hybrid) Fock 
components in DIS is discussed and the gluon contribution to structure func- 
tions and proton spin is emphasized. The last section concludes the general 
picture of the QCD dynamics with the discussion of the selfconsistent calcu- 
lation of the FC. 
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2 The QCD vacuum structure. Stochastic 
vs coherent. 



The basic quantity which defines the vacuum structure in QCD is the field 
correlator (FC) 

D ( - n \x 1 , ...x n ) = (Fp in (x 1 )$(x ll X2)F^ n (x 2 )$(x2 1 x 3 )...F l ^ Vn (x n )$(x ni x 1 )) 1 

(1) 

$(x,y) = Pexpig / A^dz^. 
Jy 

The set of FC (JTJ) for n = 2,3, ... gives a detailed characteristic of vacuum 
structure, including field condensates (for coinciding x\ = x 2 = ■■■x n ). On 
general grounds one can distinguish two opposite situations: 1) stochastic 
vacuum 2) coherent vacuum. In the first case FC form a hierarchy with the 
dominant lowest term (x±, x 2 ) (x\ — x 2 ), while higher FC are fast 

decreasing with n. We shall call this situation the Gaussian Stochastic Ap- 
proximation (GSA). In the second case all FC are comparable, and expansion 
of physical amplitudes as the series of FC is impractical. This is the case for 
the gas/liquid of classical solutions, e.g. of instantons, magnetic monopoles 
etc. The physical picture behind the situation of nonconverging FC series 
is that of the coherent lump(s), when all points in the lump are strongly 
correlated. 

To understand where belongs the QCD vacuum one can start with the 
Wilson loop in the representation D of the color group SU(3), 

W D (C) = (tr D exp(zg [ dz^T™)) (2) 

The Stokes theorem and the cluster expansion identity allow to obtain 
the basic equation, which is used in most applications of the FCM (for more 
details see 0) 

W(C) = expJ2^j- I D^ n \x 1 ,...x n )da^ Ul (x 1 )...da tlnVn (x n ). (3) 
n n\ J 

Here integration is performed over the minimal surface S m i n inside the 
contour C defined in (J2J) and D n is the so-called cumulant or the connected 
correlator, obtained from the FC Eq.Q by subtracting all disconnected av- 
erages. From (J3J) one easily obtains that the Wilson loop has the area- 
law asymptotics, W(C) ~ exp(—aS min ), for any finite number of terms 
n ma x\ n < n max in the exponent ©. 
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The string tension is expressed through D^ n \ 

o = - J D (2 \ Xl - x 2 )d 2 {x 1 - x 2 ) + 0(D {n \n > 4) = a 2 + a 4 + ... (4) 

Eq.(@J) has several consequences: 1) confinement appears naturally for 
n = 2, i.e. in the GSA 2) the lack of confinement can be due to vanishing of all 
FC, or due to the special cancellation between the cumulants, as it happens 
for the instanton vacuum jl], 3) for static quarks in the representation D of 
the color group SU(3), the string tension a 2 is proportional to the quadratic 
Casimir factor ( the Casimir scaling) 

( D ) {fund) ^(2) L / 2 , i 2 , o , o \ / K \ 

a 2 = (fund) ** > C D =o(^ + V +3/X + 3^). (5) 

C D 

However for larger n, n > 4 the Casimir scaling is violated: 

aP = cuCff + a 2 {C^f + a 3 C§ ] + ... (6) 

It is remarkable that perturbative interaction of static quarks (r) satisfies 
the Casimir scaling to the order 0(g 6 ) considered so far so the total 
potential V^ D \r) = V,^l{r) + a^r + const is also Casimir scaling, if GSA 
works well. 

This picture was tested recently on the lattice jE] and confirmed the 
Casimir scaling with the accuracy around 1% in the range 0.1 < r < 1.1. fm. 
The full theoretical understanding of this fundamental fact is still lacking, 
both for the perturbative part and for the string tension. On the pedestrian 
level the Casimir scaling and the quadratic (Gaussian) correlator dominance 
implies that the vacuum is highly stochastic and any quasiclassical objects, 
like instantons, are strongly suppressed in the real QCD vacuum. The vac- 
uum consists of small white dipoles of the size T g made of neighboring field 
strength operators. The smallness of T g might be an explanation for the 
Gaussian dominance since higher correlator terms in a are proportional to 
(FT 2 ) n (T 2 )-\ where F is the estimate of the average nonperturbative vac- 
uum field, F ~ 500 (MeV) 2 . 

Lattice calculations of FC have been done repeatedly during last decades, 
using cooling technic [7] and with less accuracy without cooling jH] . (Recently 

lr rhe diagrams of the order 0(g s ) violating the Casimir scaling have been found by 
W.Wetzel (to be published) 
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another approach based on the so-called gluelump states was exploited on 
the lattice [§] and analytically ^U], which has a direct connection to FC). 

The basic result of |7] is that FC consists of perturbative part 0(1/ x 4 ) at 
small distances and nonperturbative part 0(exp(— x/T g )) at larger distances 
with T g in the range T g = 0.2 fm (quenched vacuum) and T g = 0.3 fm ( 2 
flavours) . 

Calculations in [3] and [HHHI], as well as sum rule estimates [TT] yield a 
smaller value, T g « 0.13 fm to 0.17 fm. This enables us in what follows to 
take the limit T g — > keeping a = const « 0.18 GeV 2 , and consider aTg as 
a small parameter of expansion, aTg <C 1. For example, the contribution of 
higher correlators in (J1J is proportional to a(aTg)^~ l ,n = 4, 6, 8, ... 



3 Hamiltonian for valence components 

There are two possible approaches to incorporating nonperturbative field 
correlators in the quark-antiquark (or 3q) dynamics. The first has to deal 
with the effective nonlocal quark Lagrangian containing field correlators ^3] • 
From this one obtains first-order Dirac-type integro-differential equations for 
heavy-light mesons light mesons and baryons ^H]- These equations 

contain the effect of chiral symmetry breaking jT21 EI] which is directly con- 
nected to confinement. 

The second approach is based on the effective Hamiltonian for any gauge- 
invariant quark-gluon system. In the limit T g — > this Hamiltonian is simple 
and local, and in most cases when spin interaction can be considered as a 
perturbation one obtains results for the spectra in an analytic form, which 
is transparent. 

For this reason we choose below the second, Hamiltonian approach ^SJ 
IT7] . We start with the exact Fock-Feynman-Schwinger Representation for 
the qq Green's function (for a review see |18j). taking for simplicity nonzero 
flavor case ^ ^ 

G^ v) = j o dsi J o ds 2 (Dz) xy (Dz) xy e- K '- K " 

(trT in ( mi - ^WriQT^imi - D 2 )) A (7) 

where Ki = J Sl driirrii + \{z^) 2 ), T in>out = 1,75, ... are meson vertices, and 
W a (C) is the Wilson loop with spin insertions, taken along the contour C 
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formed by paths (Dz) xy and (Dz) xy , 

W a (C) = P F P A exp(ig f A^dz^x 

J c 

xexp(gj^ a^JF^dn-g J a^F^d^). (8) 

The last factor in (j§J) defines the spin interaction of quark and antiquark. 
The average {W )a in © can be computed exactly through field correlators 
(F {!)... F{n)) a, and keeping only the lowest one,(F(l)F(2)), which yields 
according to lattice calculation accuracy around 1% 5,, one obtains 

(W a (C)) A ~eM-hf ds^(l) [ ds Xa (2)+ 

+ E H a^dTjKF^m^))). (9) 

iij=1 Jo " Jo 

The Gaussian correlator (F MI/ (l)Fx CT (2)) = D^\ a (l,2) can be rewritten 
identically in terms of two scalar functions D(x) and D\(x) [2], which have 
been computed on the lattice [7j to have the exponential form D(x), Di(x) ~ 
exp(— \x\/T g with the gluon correlation length T g « 0.2 fm. 

As the next step one introduces the einbein variables fii and u; the first 
one to transform the proper times Sj, Ti into the actual (Euclidean) times 
ti = z| . One has [T7j 

2^) = ^ tfeipSW)^ = const y Dynit^z^y. (10) 

The variable enters in the Gaussian representation of the Nambu-Goto 
form for SVnm and its stationary value vq has the physical meaning of the 
energy density along the string. In case of several strings, as in the baryon 
case or the hybrid case, each piece of string has its own parameter i/W. 

To get rid of the path integration in Q one can go over to the effective 
Hamiltonian using the identity 

G q g(x,y) = (x\exp(-HT)\y) (11) 

where T is the evolution parameter corresponding to the hypersurface chosen 
for the Hamiltonian: it is the hyperplane z^ = const in the cm. case |17j . 
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The final form of the cm. Hamiltonian (apart from the spin and pertur- 
bative terms to be discussed later) for the qq case is [T71 [TH] 



2 Jo u((3) Jo 2 H K ' 

Here ( = (pi + J flvdfi) / + ^2 + So fivdfl) and and v{0) are to be 
found from the stationary point of the Hamiltonian 

9H °\ -n dH °\ -n hi) 

a^ 1 ^" ' ^7>— (0) -°- (13) 

Note that H contains as input only m^mj and a, where are current 
masses defined at the scale of 1 GeV. The further analysis is simplified by the 
observation that for L = one finds z/°) = ar from (j!3|) and /i, = \/m 2 + p 2 , 
hence ifo becomes the usual Relativistic Quark Model (RQM) Hamiltonian 



h (l = 0) = \K 2 + p 2 + or - ( 14 ) 



i=i 



But is not the whole story, one should take into account 3 additional 
terms: spin terms in (jHJ) which produce two types of contributions: self-energy 
correction HC 



H self = ±^, A^ = -%K), ,(0) = 2(1 + ^#)^ 1*0.9 

(15) 

and spin-dependent interaction between quark and antiquark H spin [21] which 
is entirely described by the field correlators D(x), Di(x), including also the 
one-gluon exchange part present in Di(x). 

Finally one should take into account gluon exchange contributions, which 
can be divided into the Coulomb part Hc ou i = — §^7^, and H ra d including 
space-like gluon exchanges and perturbative self-energy corrections (we shall 
systematically omit these corrections since they are small for light quarks 
to be discussed below). In addition there are gluon contributions which are 
nondiagonal in number of gluons n g and quarks (till now only the sector 
n g = was considered) and therefore mixing meson states with hybrids and 
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glueballs [221; we call these terms H mix and refer the reader to [22| an d the 
cited there references for more discussion. Assembling all terms together one 
has the following total Hamiltonian in the limit of large N c and small T g (for 
more discussion see j2H]) : 

H = H + H se if + H spin + HqouI + H ra d + H mix . (16) 

We start with H = Hr + H string . The eigenvalues M of Hr can be given 
with 1% accuracy by [21j 

Ml «8<rL + 47r<r(n + ^) (17) 

where n is the radial quantum number, n = 0,1,2,... Remarkably M ~ 
4/i , and for L = n = one has /x (0,0) = 0.35 GeV for a = 0.18 GeV 2 , 
and /j,Q is fast increasing with growing n and L. This fact partly explains 
that spin interactions become unimportant beyond L = 0, 1, 2 since they are 
proportional to dridr 2 ~ 4 * cfcicfo (see (J0J) and |2Sj)- Thus constituent 
mass (which is actually "constituent energy") /i is "running". The validity 
of socially accepted " constituent mass" is confirmed by its numerical 

value given above, the spin splittings of light and heavy mesons [2H] and 
by baryon magnetic moments expressed directly through fio, and being in 
agreement with experimental values 



4 Hamiltonian and bound states of valence 
gluons 

We now come to the gluon-containing systems, hybrids and glueballs. Refer- 
ring the reader to the original papers (23-1211 one can recapitulate the main 
results for the spectrum. In both cases the total Hamiltonian has the same 
form as in (|16|). however the contribution of corrections differs. 

For glueballs it was argued in [29J that Hq fll2j) has the same form, but 
with m,; = and a — > a a dj = \o while H se if = due to gauge invariance. 

Thus one can retain in (fTo^) only two main terms: H = Hq + H sp i n 
while Hcoui was argued to be strongly decreased by loop corrections. The 
calculation in [21] was done for two-gluon and three-gluon glueball states and 
results are in surprisingly good agreement with lattice data for both systems 
(no fitting parameters have been used in [29 J. 
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We now coming to the next topic of this talk: hybrids and their role in 
hadron dynamics. We start with the hybrid Hamiltonian and spectrum. This 
topic in the framework of FCM was considered in j23 I2E] The Hamiltonian 
H for hybrid looks like HH |2E| 

= ^1 + ^1 + a + ^+Mg + P|±P^ + a ±\r g - r,| + H str . (18) 
2fn 2/x 2 2 2/i £f 3 

Here p^, p^ are Jacobi momenta of the 3-body system, H se if is the same 
as for meson, while H spin and H Cou i have different structure [2%] . 

The main feature of the present approach based on the BPTh, is that 
valence gluon in the hybrid is situated at some arbitrary point on the string 
connecting quark and antiquark, and the gluon creates a kink on the string 
so that two pieces of the string move independently (however connected at 
the point of gluon). This differs strongly from the flux-tube model where 
hybrid is associated with the string excitation as a whole, but has a strong 
similarity to the treatment of gluons in the framework of the Lund model 

Results for light and heavy exotic 1 + hybrids are also given in [23j and 
are in agreement with lattice calculations. Typically an additional gluon in 
the exotic (L = 1) state "weights" 1. 2-=-1.5 GeV for light to heavy quarks, 
while nonexotic gluon {L = 0) brings about 1 GeV to the mass of the total 
qqg system. Let us now consider the hybrid spectrum in more detail. First 
of all we use for that 3-body problem the hyperspherical method, which 
works with accuracy of few percent (311 E21- Then the whole spectrum is 
classified by the grand angular momentum K = 0, 1, 2, which is actually 
an arithmetic sum of all partial pair angular momenta in the system qqg. 
The lowest K — states can be formed from the s-wave qq pair and s-wave 
valence gluon g, which gives the p + g and 7r + g systems, and vectors imply 
spin-one particle. In this way one obtains the classification 

K = 0, (tt + g) - 1+- (p + g) - (2++, 1++, 0++) 

K = 1, (tt + (V x g)) - 1-, (p + (V x g)) - (2-+, r+, 0-+). 

The eigenvalues of the Hamiltonian (|18j) are easily obtained for the light 
quarks using the hypercentral (lowest K) approximation - 

M(K = 0) = 1.872 GeV 
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M(K = 1) = 2.45 GeV 

M(K = 2) = 2.90 GeV 

M(K = 3) = 3.27 GeV. 

Here the (negative) self-energy part of quarks H se if is already added to 
masses. One has also in addition the color Coulomb part Hcoui and spin- 
dependent part H spin , which contribute approximately (H Co ui) ~ —0.2 GeV 

( ~ 2 \ 

and for spin-spin interaction approximately 0.08 GeV —1 , where num- 

bers inside brackets refer to 1 + 1 = 0, 1, 2 respectively from top to bottom. 

As a result, neglecting rather small contribution from H ra ^ and H mix in 
(JUIJ), and the string correction, taking into account the moment of inertia of 
the string [T7] . yielding around (-100 MeV), one has the approximate mass 
estimates for the lowest mass states of K multiplets: 

M low (K = 0) = 1.42 GeV 

M low (K = 1) = 1.9 GeV 
M low (K = 2) = 2.45 GeV. 

The mass value Mi 0W (K = 1) agrees well with the lattice results for the 
mass of the 1~ + state |53j . The recent unquenched calculation yields the 
value which is somewhat lower. 

One should stress that the hybrid states, which start at the mass around 
1.4 GeV, have a high multiplicity which grows exponentially with mass, as 
well as excited string states in bossonic string theory |22J. This fact has 
a very important consequence for high-energy processes, where the hybrid 
excitation is argued to be the dominant physical mechanism. 



5 Hamiltonian and Fock states 

As was mentioned above the QCD Hamiltonian is introduced in correspon- 
dence with the chosen hypersurface, which defines internal coordinates {^} 
lying inside the hypersurface, and the evolution parameter, perpendicular to 
it. Two extreme choices are frequently used, 1) the cm. coordinate system 
with the hypersurface £4 = const., which implies that all hadron constituents 
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have the same (Euclidean) time coordinates x 4 = const, i = 1, ...n, 2) the 
light-cone coordinate system, where the role of £4 and is played by the 
x + ,x^ components, x + = x °^' 9 ■ 

To describe the structure of the Hamiltonian in general terms we first as- 
sume that the bound valence states exist for mesons, glueballs and baryons 
consisting of minimal number of constituents. To form the Fock tower of 
states starting with the given valence state, one can add gluons and qq pairs 
keeping the J PC assignment intact. At this point we make the basic sim- 
plifying approximation assuming that the number of colors N c is tending to 
infinity, so that one can do for any physical quantity an expansion in powers 
of 1/N C . Recent lattice data confirm a good convergence of this expansion 
for N c = 3, 4, 6 and all quantities considered IHH] (glueball mass, critical 
temperature, topological susceptibility etc.). 

Then the construction of the Fock tower is greatly simplified since any 
additional qq pair enters with the coefficient 1/N C and any additional white 
(e.g. glueball) component brings in the coefficient 1/iV^. In view of this in 
the leading order of 1/N C the Fock tower is formed by only creating additional 
gluons in the system, i.e. by the hybrid excitation of the original (valence) 
system. Thus all Fock tower consists of the valence component and its hybrid 
equivalents and each line of this tower is characterized by the number n of 
added gluons. Then, the internal coordinates {£}„ describe coordinates and 
polarizations of n gluons in addition to those of valence constituents. 

We turn now to the Hamiltonian H, assuming it to be either the total 
QCD Hamiltonian Hqcd, or the effective Hamiltonian H^ e ^\ obtained from 
Hqcd by integrating out short-range degrees of freedom. We shall denote the 
diagonal elements of H, describing the dynamics of the n-th hybrid excita- 
tion of s-th valence state (s = m{f f}, gg,3g,b{f 1/2/3} for mesons, 2-gluon 
and 3-gluon glueballs and baryons respectively with fa denoting flavour of 
quarks) as H^. For nondiagonal elements we confine ourselves to the lowest 
order operators H^ n+1 and H^_ ln describing creation or annihilation of one 
additional gluon, viz. 




(19) 



H g2g = f r 6c /(ax - d ua y) a yj\ 



(20) 



and we disregard for simplicity the terms H g3g . 
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As it is clear from (|19|). (|2(Jj). the first operator refers to the gluon creation 
from the quark line, while the second refers to the creation of 2 gluons from 
the gluon line. In what follows we shall be mostly interested in the first 
operator, which yields dominant contribution at large energies, and physically 
describes addition of one last cross-piece to the ladder of gluon exchanges 
between quark lines, while (j2Ti|) corresponds in the same ladder to the a s 
renormalization graphs, and to the graphs with creation of additional gluon 
line. 

The effective Hamiltonian in the one-hadron sector can be written as 
follows 

H = H<® + V (21) 
where is the diagonal matrix of operators, 

HW={H<g,H8,H<$,...} (22) 

while V is the sum of operators ()19|) and (|20p. creating and annihilating 
one gluon. In (}2*2*|) is the Hamiltonian operator for what we call the 
"n-hybrid", i.e. a bound state of the system, consisting of n gluons together 
with the particles of the valence component. In this way the n-hybrid for the 
valence p-meson is the system consisting of qq plus n gluons "sitting" on the 
string connecting q and q. 

Before applying the stationary perturbation theory in V to the Hamil- 
tonian (|21|). one should have in mind that there are two types of excita- 
tions of the ground state valence Fock component: 1) Each of the operators 
Hnni n = 0, 1, ... has infinite amount of excited states, when radial or orbital 
motion of any degree of freedom is excited, 2) in addition one can add a gluon, 
which means exciting the string and this excitations due to the operator V 
transforms the n — th Fock component ijj^ into tpn+i- 

The wave equation for the Fock tower ty N {P, £} has the standard form 

H^ N = (H {0) + V)V N = E N q N , (23) 
or in the integral form 

y N = *§ ) -G®VV N (24) 
where is diagonal in Fock components, 
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and iffff is the eigenfunction of H^°\ 

H^P = (26) 

and since is diagonal, has only one Fock component, = ip n {P, 
n = 0,l,2,..., and the eigenvalues E$ contain all possible excitation energies 
of the n-hybrid, with the number n of gluons in the system fixed, 

E$ ) =EW{P) = y /P* + MZ w . (27) 

Here {k} denotes the set of quantum numbers of the excited n-hybrid. 

From ()24j) one obtains in the standard way corrections to the eigenvalues 
and eigenfunctions. 

As a first step one should specify the unperturbed functions ' , intro- 
ducing the set of quantum numbers {k} defining the excited hybrid state for 
each n-hybrid Fock component ?/> n (P{£}); we shall denote therefore: 

=VW P >{£}n), 71 = 0,1,2,... (28) 

The set of functions ip n {k} with all possible n and {k} is a complete set 
to be used in the expansion of the exact wave-function (Fock tower) ty N : 

= E C rn{k}^m{k}- (29) 
m{k} 

Using the orthonormality condition 

/ ^J^W^ = 5 rnJ{k}{p} (30) 

where dT implies integration over all internal coordinates and summing over 
all indices, one obtains from (JT5|l an equation for c m ^} and E N , 

Cn{p}( E N - E n{p}) = C m{k} V n{p},m{k} ( 31 ) 
m{k} 

where we have defined 

V n { p }, m {k} = J ip^ {p} Vi> m {k}dT. (32) 

Consider now the Fock tower built on the valence component ip v i K \, where 
v can be any integer. For = 0{0} this valence component corresponds 
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to the unperturbed hadron with minimal number of valence particles. For 
higher values of the Fock component ip v i K \ corresponds to the hybrid 

with v gluons which after taking into account the interaction is "dressed 
up" and acquires all other Fock components, so that the number N in (129)1 
contains the "bare number" as its part iV = ... (at least for small 

perturbation V"). 

One can impose on \I/jv the orthonormality conclusion 

J ^ M dT = J2 = <W (33) 

m{k} 

Expanding now in powers of V, one has 

N(u{*}) _ r r , Ml) N{2) 

C m{k} ~ OrrwO {k}{K} + C m{fc} + C m{fc} + ... (34) 

E N ( U { K }) = E^ k} + E$ + E$ + ... (35) 

It is easy to see that E$ = 0, while for one obtains from (|3*T|) the 
standard expression 

JV(1) _ Vn{p},u{K} 
°n{p} - p (0) p (0) • W 

In what follows we shall be interested in the high Fock components, u+l, {k}, 
obtained by adding I gluons to the valence component z/{k}. Using ()31J) and 
one obtains 



c 



N(u{k}) _ \p Vy+ljky^+l-ljk!} K+l-l{fci} ,v+l-2{k 2 } 

u+i,{k} ~ 2^ Tjo) 7f(o) T^o) ' 



V, 



v+l{k,},v{n} n( ^ / i + 2 



P (0) p(0) 



+ 0(V t+2 ). (37) 



Since V is proportional to g, one obtains in (}3"4"jl the perturbation series 
in powers of a s for c N and hence for (|29|1 . One should note that a s (Q 2 ) 
is the background coupling constant, having the property of saturation for 
positive Q 2 j3] and the background perturbation series has no Landau ghost 
pole and is defined in all Euclidean region of Q 2 . 

The estimate of the mixing between meson and hybrid was done ear- 
lier in the framework of the potential model for the meson in ^7\. In 
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the mixing between hybrid, meson and glueball states was calculated in the 
framework of the present formalism and we shortly summarize the results. 
One must estimate the matrix element (}3*2"|) between meson and hybrid wave 
functions taking the operator V in the form of ([19)1 . where the operator of 
gluon emission at the point (x, 0) can be approximated as 

1 



x [exp(zk • x - ^t)ej t A) c A (k) + ej, A) c£(k) exp(-zk ■ x + fat)] (38) 

Omitting for simplicity all polarization vectors and spin-coupling coeffi- 
cients which are of the order of unity, one has the matrix element 

V Mh = -J= J <p M (T)^+(0, r)d 3 r (39) 

where ipm{ t ), M VVi( r i> r 2) are meson and hybrid wave functions respectively, 
and in (|39|) it is taken into account that the gluon is emitted (absorbed) from 
the quark position. 

Using realistic Gaussian approximation for the wave functions in (J39)) one 
obtains the estimate 



V Mh ^9- 0.08 GeV. (40) 

A similar estimate is obtained in |22] for the hybrid-glueball mixing matrix 
element, while the meson-glueball mixing is second-order in (J4(Jj) . 
Hence the hybrid admixture coefficient (|36|) for the meson is 

ri VMh _ VMh , . 1 n. 

and for the ground state low-lying mesons when AMm/j ~ 1 GeV it is small, 
Cjsth ~ 0.1 — 0.15, yielding a 1-2% probability. However for higher states in 
the region M M > 1.5 GeV, the mass difference AM Mh of mesons and hybrids 
with the same quantum numbers can be around 200 MeV, and the mixing 
becomes extremely important, also for meson-glueball mixing, which can be 
written as 

n sr VwhVhG Ur) s 

° MG ~ \ AM Mh AM hG (42) 

and Vj^h ~ VhG- It is clear that the iterative scheme described above can be 
useful only because hybrid excitation by one additional gluon "costs" around 
1 GeV increase in mass, hence the coefficient fl3fij) can be small. 
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6 Hybrid states and DIS 

As was stressed in the previous section, in the large N c limit the higher Fock 
components which are excited by the external current (or incident hadron) 
are the multihybrid (or n- hybrid), states. It is convenient to consider these 
states in the light-cone formalism, following the line of derivation and most 
notations in [3*T] . 

Consider the n-hybrid with quark at the point zft , antiquark at the point 
zjp and gluons at the points z^ =1 n . We also define = —z( l ~ 1 \ with 

z (0) = z (a) and z (n+l) = z (b)_ 

The action is 

A = K + aS min (43) 
where the kinetic operator K and the minimal area S m i n can be written as 

K=^ + ^ + \l T dz + [p a ((z{ a) ) 2 + 2^ a) ) + Ho((z?) 2 + 2i W ) 



2p a 2/i fc 2 Jo 



+ 5>((4 i) ) 2 + 2ii i) )] (44) 



i=l 

\ pT n +l pi 



oS min = - cfe+ V / df3i 
2, Jo t^i Jo 



V; 



(45) 



( W 'tt)2 J 

and we have defined 

™f = ^r i] (l - Pi) + &z®, w {t) = z^\l - pi) + (46) 
u ,'W = z (0_ iZ «-i)= p (i). (47) 



As in [37] we introduce the total momentum P + , 

n+l -i 
10 



n n+l -i 

P+ = Y,Vi + J2 / V * d P + + fib- (48) 

■ 1 ■ i <j o 



i=l i=l 

At this point one can make a new important step and introduce the parton's 
quota Xi of the total momentum P + to be associated with the Feynman 
variables X{, i — 1, ...,n, x a and £&, which can be written as 

_ tM + Jo ViW + Jo ~ , AQ x 
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X a = Xq = (fl a + [ - ff)d(3)/P. 



Jo 



(50) 



X\, — X n -\-\ 



+ f u n+1 (3df3)/P + 
Jo 



(51) 



One can notice that Xi consists of three pieces: 1) the (+) -component of 
momentum of the valence gluon (/ij), 2) the (+) momentum of the preceeding 
piece (i) of string weighted with the factor f3 which takes into account that 
the string (i) is deformed by the motion of the gluon (i) while another end of 
the string is fixed 3) the (+) momentum of the string (i + 1) weighted with 
the factor (1 — 0) taking into account motion of the (i ^ 1) string due to the 
z-th gluon. Note that the parameter (3 in all strings (i), i — 1, ...n+ 1 grows 
in one direction, e.g. from the left to the right. 

In this way the momentum of each piece of the string (i) is shared by 
two adjacent gluons: (i — 1) and (z), so that each parton quota Xi contains 
momentum of the parton (gluon or quark or antiquark) itself and of pieces 
of adjacent strings. 

These results imply several nontrivial consequences. First of all, one can 
see that the einbein factors /!,, which played in the cm. system the role 
of constituent mass (energy) of gluon and quarks, and played the role of 
energy density along the string, in the I.e. system they enter directly the 
Feynman parameters of gluons. In this way one can for the first time see 
the connection of the standard constituent quark (gluon) picture with the 
parton picture and calculate as in (}50|) . (j5Tf parton parameters through the 
(Lorentz boosted) constituent energies of quarks and gluons. Secondly in the 
I.e. wave-function of the n-hybrid the average values of \ii and Vi are equal for 
large n, Jii = v^i = 1, 2, ...n, while for n = 1 one obtains jj, g = y/2p, g = \f2jiq 



Hence gluons carry more momentum on average than quarks in the n- 
hybrid state. Therefore one expects that in DIS at large enough energy 
when the n-hybrid component of the hadron wave-function is excited, the 
contribution of gluons to the momentum sum rule and to the proton spin 
should be dominant. This expectation is consistent with the experimental 
data. Thus in the neutrino-isoscalar scattering the momentum sum rules for 
the quark part of F 2 (x, Q 2 ) yield jHO] 0.44± 0.003, which implies that gluons 
carry more than 50% of the total momentum. 



[38J. 



17 



For the proton spin one has the relation [3U] 

J = iAEM + L q Qi) + J a {fi) (52) 

where the quark sigma-term experimentally is AE(/x = 1 GeV) = 0.2 ± 0.1, 
and the most part of the difference between the l.h.s. and the r.h.s. is 
presumably due to the gluon spin contribution J g {n). The detailed estimates 
of hybrid contributions to DIS and high-energy scattering will be published 
elsewhere |3*8] . 

7 Conclusions 

It was explained above that the Field Correlator Method is a powerful tool 
for investigation of all nonperturbative effects in QCD. In particular it pro- 
vides a natural mechanism of confinement, compatible with all lattice data, 
and explains the close connection of confinement and chiral symmetry break- 
ing. The spectrum of mesons, glueballs and hybrids is calculated with a, a s 
and current masses as the only fixed parameters used and this spectrum 
is in good agreement with lattice data and experiment. The latest de- 
velopment concerns the dominant role of hybrids in DIS and high-energy 
scattering and here the first qualitative results are consistent with exper- 
imental evidence. The author is grateful to the organizers of the Pomer- 
anchuk International Conference for their excellent job, and to A.M.Badalian, 
K.G.Boreskov, A.B.Kaidalov and O.V.Kancheli for many stimulating discus- 
sions. 

The partial support of the INTAS grants 00-110 and 00-366 is gratefully 
acknowledged. 
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